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I present results of recent work for the UKQCD collaboration on the light hadron spectrum using a non- 
perturbatively 0(a)-improved Wilson action with two degenerate flavours of dynamical quarks on a 16 3 x 32 
lattice. Values of the bare gauge coupling and bare dynamical quark mass were chosen to keep the lattice spacing 
fixed at two distinct values in physical units as determined through the Sommer scale parameter, ro. I also include 
results for a quenched ensemble with a lattice spacing matching one of these dynamical sets and results from an 
exploratory run at lighter bare quark mass. 



1. Introduction 

The prediction of the spectrum of light hadron 
masses is one of the fundamental aims of lattice 
QCD. This has already been done extensively in 
the quenched approximation by several collabo- 
rations. The results for the quenched spectrum 
have been found to lie within 10% of the experi- 
mental values. This is encouraging and motivates 
a full calculation with dynamical quarks. The re- 
sults of such a study would test the appropriate- 
ness of QCD as the theory of the strong interac- 
tion. 

In this paper I present work on the light hadron 
spectrum for the UKQCD Collaboration using a 
non-perturbatively 0(a)-improved Wilson action 
with two degenerate dynamical quarks. Much of 
this work can be found in a recent paper on un- 
quenching effects by the UKQCD Collaboration 
The three dynamical simulations and one 
quenched simulation used to study unquenching 
were performed at values of f3 and K sea (where ap- 
propriate) falling on a trajectory of fixed lattice 
spacing. This was done to minimise variations 
due to residual discretisation errors and finite vol- 
ume effects. In addition to these four simulations 
I include results from two dynamical simulations 
lying on a different fixed lattice spacing trajectory 
and one exploratory run at lighter quark mass to 
test the limits of the current algorithms. The 



fixed lattice spacing trajectories were determined 
using the matching techniques in || and the scale 
set using ro in all simulations. 

2. Simulation Details 

All of the simulations discussed here were per- 
formed on a 16 3 x 32 lattice. For the dynam- 
ical simulations we used the Wilson gauge field 
action together with the Sheikholeslami-Wohlert 
0(a)-improved Wilson gauge- fermion action ||. 
The gauge configurations were generated using 
the Hybrid Monte Carlo algorithm Q. 

Quark propagators were calculated using the 
0(a)-improved gauge- fermion term for both dy- 
namical and quenched simulations. We used the 
final published values of the 0(a)-improvcmcnt 
coefficient c sw , as determined by the Alpha col- 
laboration in the generation of dynamical 
gauge fields and quark propagators. 

In what follows I distinguish the values of the 
hopping parameter, k, used in the generation 
of gauge configurations and quark propagators 
by K S ea and n va i respectively. The cases where 
K-vai 7^ Ksea were studied to allow a partially 
quenched analysis to be performed. 

I summarise the simulation parameters in Ta- 
ble 1. The value of L/tq is included as an indi- 
cator of the onset of finite size effects. The pa- 
rameter, ro, was determined from calculations of 
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Table 1 



Summary of simulation details, statistics and specific observables with statistical errors only. Gauge 
configuration for the (5.2,0.13565) dataset have been separated by 20 trajectory steps ( as opposed to 40 
for the other datasets) since statistics were low. 
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the inter-quark potential and was used to set the 
scale. I also include MpgjMy to indicate how 
light the simulations are. This ratio was calcu- 
lated at K S ea — Kval • 

It can be seen that all the simulations were per- 
formed with relatively coarse lattices. This was 
to avoid the appearance of finite size effects and 
was limited by the available computer resources. 
Previous analysis by UKQCD || has shown that 
finite size effects in the light hadron spectrum are 
negligible provided 



lators at time, t, 

N 

C(t) = ]T A n ( e - M » 4 + e - M "( T - 4 >) . (2) 
n=0 

Here T is the temporal extent of the lattice, A n is 
a time independent factor and M n is the mass of 
the n th excited state with the ground state(n = 0) 
being the mass we wish to calculate. For baryons 
a similar ansatz was used but with the second 
exponential term omitted. 



L/r > 3.2 



(1) 



3. Fitting Hadron Correlators 

In this work I have concentrated on determin- 
ing the masses of the pseudoscalar and vector 
mesons and the nucleon and delta baryons. The 
hadronic interpolating operators used in the cal- 
culation of the necessary two-point correlation 
functions were cither local or smeared using the 
fuzzing procedure of [7]. I use a nomenclature 
where an FL correlator is one that is fuzzed at 
sink and local at source. 

It should be noted that hadron correlators were 
calculated for all the possible combinations K va i 
as given in Table |l| although baryons with non- 
degenerate constituent quark masses have not 
been analysed yet. 

The meson masses were determined from their 
respective correlation functions by correlated \ 2 
minimisation with the following ansatz for corre- 



o LL 

□ FL 

O FF 

A LL.FL.FF 



J L 

10 



Figure 1. Vector fit masses for the (3 = 
5.25, K sea = 0.1352 dataset with n va i — n sea ver- 
sus the time range minimum with t max fixed at 
15. 
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Table 2 

Preliminary hadron mass results in lattice units for the (5.2,0.13565) and (5.25,0.1352) datasets with 
statistical errors only. The values for n va i^ are only relevant for baryonic results. 
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In addition to the fits to single correlators, si- 
multaneous fits were performed to LL, FL and FF 
correlators. Not only did this allow fitting to a 
common mass but the amplitude factors, A n , for 
the three correlators could be written in terms of 
two common factors, further constraining the fit 
and reducing the total number of parameters. 

The hadron mass analysis proceeded in all cases 
by performing fits over a series of time ranges. As 
an example Figure |l| shows the dependence of the 
fitted ground state mass for a vector meson on 
the minimum time of the fit range, t m i n . Final 
fit ranges were chosen when the goodness of fit 
parameter, Q, indicated the ansatz was appropri- 
ate and when some stability in the fit parame- 
ters with respect to small changes in the fit range 
had been reached. The statistical errors on the 
fit parameters were calculated by the bootstrap 
method. 



4. Results 

The fit results for all but the (j3 = 5.25, K sea = 
0.1352) and (/3 = 5.2, n sea = 0.13565) datasets 
can be found in (lj. I present preliminary re- 
sults for these two remaining datasets in Table 
H Fits to mesonic correlators were performed for 
all combinations of the valence quark masses, as 
given by K va ii and K va n- Baryonic fits have only 
been performed with degenerate valence quarks. 

All of these results are from simultaneous fits 



to LL, FL and FF correlators. The results for 
the (5.2, 0.13565) dataset and the mesonic fits for 
the (5.25,0.1352) dataset were produced with a 
ground state only fit ansatz (i.e. N — in equa- 
tion |^). The remaining baryonic fits were per- 
formed with the first excited state contribution 
included. 
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